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Abstract 

Oh! 

D . The screened quasi-relativistic potential is used for describing spin-orbit splitting 

' in ^Pj waves of quark-antiquark system. Fermi-Breit equation is solved numerically 

in configuration interaction approximation. This approximation takes into account 
^ . the mixing of wave functions up to fifth order and corrects substantially pertur- 

^ I bation calculations. The nature of potentials Lorentz transformation property is 

elucidated. Very good quantitative results for bb and cc quarkonia and quite ac- 
ceptable qualitative characteristics as well as for systems with unequal masses are 
obtained for uu. 
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Today it seems evident that quark potential model gives a rather good description of 
spin-averaged mass spectrum of hadrons, considered as a system of quarks (see e.g. |l|] 
and literature cited there). Usually non-relativistic Cornell as well as oscillator potential 
with Coulumb-like one-gluon exchange or other power-law confinement terms is used. In 
this work we try to extend this approach to incorporate the second order spin-terms in 
two-quark Fermi-Breit equation for evaluating the spin-orbit splitting. Instead of calcu- 
lating as usual spin-terms in first-order perturbation approximation the expansion of the 
total wave function into a basic set of non-perturbative solutions up to the fifth order of 
configurationally interacting states is carried out. 

The influence of a specific form of potential on a spin-effects (like Cornel, oscillator, 
screened etc.) which generally give the same results for spin-averaged data, the way of 
accounting of both long-range non-perturbative and short-range (perturbative) parts in 
these calculations, the Lorentz structure of different parts of potentials all these questions 
remain to be unsolved yet. 

Our main goal is to clarify some aspects of these questions in framework of config- 
uration interaction (CI) approach [H. This method does not need the assumption that 
the coupling constant is to be small, the assumption which is required by perturbation 
method. Since in quark potential case it is not so the use of perturbation method looks 
very dubious. 

Let us suggest that the static quark potential has vector and scalar property of Lorentz 
transform: 

Vnr (r) = Vv (r) + Vs (r) . (1) 

Index "v" means that potential is a 4-th component of operator p^, index "s" means 
that the potential is scalar. 

Following many authors we assume the admixture of vector-scalar screened potential 

Vv ir) = VoGE + e ■ Vconf, (2) 
Vs{r) = {l-e)-Vconf, (3) 

where 

Vqge— is the one gluon exchange potential, 
Vconf— is the confinement part of the potential. 

where e is mixing constant. Here the Lorentz nature of the one-gluon and the con- 
fining potential is different, the one-gluon potential being totally vector while the con- 
fining potential is a vector-scalar mixture. This choice seems to be reasonable since 
non-perturbative vertex corrections are important at small i.e. large distances and 
naturally are coupled only with the long-range term in potential. Very interesting rewiew 
was done in the work of Nora Brambilla conserning the choice of the interaction potential 
in @. 

In our case the choice of potential itself is dictated by the consideration of most ac- 
curate description of averaged mass spectrum, and here Chikovani-Jenkovszky-Paccanoni 
(CJP) type potential 0. 

Our approach is based on the model of a non-perturbative gluon propagator, which 
was recently proposed by Chikovani, Jenkovszky and Paccanoni (CJP) [Q. In this model 
the propagator has the form 
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and the corresponding potential assumes the form |^ 

N fi'^ /-, r i\ IStt exp [-Mr] , , 

V^(r) = ^^(1-exp -H) 7^ — ^- (5) 

^ ^ 67r/i^ ^ 25 r-ln(6 + l/(Arf) ^ ^ 

In fact for accelerating numerical calculations we use simple form ^ for one-gluon- 
exchange type term with QCD-like of asymptotic freedom 

/ N 127r 1 

"^"^-^^^ lnM/(Ar) 

where A is taken to be equal to A = 0.14 GeV . The initial value of was defined via 
calculating qq— masses. 

At the finally, we present scalar and vector parts of the potential (0,^ in form 

W(r) = -^, (7) 
r 

Recently screened potential of type (|^,|^) was successfully applied to a description 
of spin-averaged meson and baryon mass-spectra P]. Recently, the study of the quark 
confining potential from QCD has been well developed. For example, the consideration 
of the gluon condense may give the confining potential a theoretical description. By 
considering the loop-diagrams of fermion in the lattice gauge calculation, E. Laermann et 
al showed that the quark confinement potential goes lower than the linear confinement 
potential, when r, the distance between quarks, becomes larger. It refiects the screening 
between valence quarks and quark sea enhances, consequently, it suppresses the strength 
of color confinement between valance quarks. According to quantum mechanics, this 
feature may reduce the energies of highly excited baryons and mesons and help us to 
understand the experimental data. 

Timo Thonhauser approached the problem of the screened potential from opposite 
site, namely of Few-Body approach to baryon structure. Few-Body potential for baryons 
usually has confinement ^ part plus screened part. Very detailed description of nucleon- 
nucleon interaction is obtained within this approach. Thonhauser applied this potential to 
meson structure and obtained quite reasonable values of averaged meson masses. Unfortu- 
nately he did not report the results concerning fine splitting. The question of comparison 
of two approaches needs further study. 

In addition, (see e. g. Gerasimov 0) the spin-orbit term has to be of short range, as is 
indicated by quantum chromodinamics (QCD), the condition which is evidently satisfied 
by CJP potential. 

In what follows we shall use the screened potential , which proved to be very good in 
describing the spin-averaged mass-spectrum of both bosons and baryons as quark systems 
0,[|l^-[jl3l and which secures the necessary fall-of the spin-dependent forces. 
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Let us start with two-body Fermi-Breit equation. We shall use nuclear system of units 

5.068 
IFm ■ 



h = c = l,lGeV = The Hamiltonian of the system has the form: 



H = Ho + W, (9) 

where 

ff„ . + ^ 4(1^) . (10) 



m is the reduced mass, 



W = Hls + Hst. (11) 



where spin-dependent potentials are given by 
• Spin-orbit interaction 



H 



LS 



1{ [((mi + ms)^ + 2mim2) ■ + (m^ - m?)^ ■ S. 
[{ml + ml) t-jt + iml-^mj)t^ Sl]jt}^ 
= ^ xp^ S+ = Si + S2, = Si - 5*2, 

t-^ = Uj{j + l)-l{l + l)-S{S + l)]. 



dVy 
dr 



Tensor terms 



(12) 



H7 



1 



S 



12 



(2«+3)(2«-l) 



12mim2 



1 dVy 

r dr 
¥2 



d-^Vv 



dr^ 



■ s 



12, 



2 ^ 



(13) 



All notations are considered to be familiar and coincide with those which were used 
by Lucha and Schoberl in their well-known monograph |jl|]. 

Lucha and Shoberl [||] has indicated that the term causes serious dimensional 

trouble because of particle falling on the center. But in a more consistent approach 
based on using Dirac equation we would have instead of previous term the following one 
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E—V+mc^ r dr 



which for r — will behave like \ and the problem is removed. In matrix 

r (ir ' T-^ 



elements, we have calculated we obtained that, there are no essential differences between 



these two results. Previously we have studied these questions in fl^ . 

The ether possibility to solve this problem is the regularisation of the potential [|l[. 



The same idea is based on the using of Vreg oc 



in the meaning of the regularized 



reg ^- r+r^ 

potential; where tq— is free parameter, which has been obtain from the comparison with 
the experimental data. So in the case of r = the problem of particle fooling is gone. 

Based on this reason we believe that our numerical results are not up set by this 
problem. 

However there are more "exact" QCD - motivated potentials, which are based on two 
loop back diagrams |T^. But by more detailed consideration it appears, that these ac- 
counts were carried out only for one-gluon exchange part. As races for systems considered 
by us, the very essential role is played by confinement. Besides use of this potential re- 
sults in increase of number of parameters, with 1 in our case up to 3. Moreover analysis 
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the obtained results with two loopback potential show that of essential improvement of 

concurrence theoretical and experimental data is not observed. 

In addition to above, indicated terms presented in Fermi-Breit Hamiltonian are the 
■ 5*2 (spin-spin) and relativistic correction terms (of p'^-order). Some authors (like [||) 

indicated that these terms are important for calculating mass spectra, other authors (like 

|1T6| ) argue that these results are purely constant. In our case we believe that for LS-mass 

differences they will not play an essential role. 
Now, we consider Fermi-Breit equation 

+ ^ ( r) = ( r*) . (14) 




There are certain difficulties in solving equation (|1^) in the case of spin-orbit-coupling. 
Using the potential ([l^) lead us to the appearance of term of order Usually authors 
use one of two possibilities: 

a) use the perturbation method; 

b) use the numerical computation. 

Both of these possibilities are inferior: first one is unacceptable because in some cases 
the "fine" splitting turns out to be not so "fine" at all, being up to 50% contribution to 
the final mass, the second one needs introducing the cut-off parameter which is highly 
undesirable. 

Here we suggest to use CI approach which was previously very successfully applied 
in atomic physics [^. The essence of this approximation is that the total wave function 
{'f) is expanded in set of eigen functions (fnoi the unperturbated Hamiltonian Hq, that 

is 

^{^) =J2(^nVn{^). (15) 

After substituting (^) into (|Tip and using eigen value we obtain the system of 
linear equations for a„ which have to be truncated for reasonably large n. 

ai {E - El - Wn) - a^Wi^ - ^sW^is - • • • - ^nW^in = 
-aiW2i + ai{E- El - W22) - 031^23 - • • • - anW2n = 

-aiWnl - a2Wn2 - asWnS - ... + aniE - E'^- Wnn) = 

where 



(16) 



W,, = {^i\W\if,). (17) 

Both the basic functions ipi and matrix elements Wij are calculated numerically. Non- 
trivial solution there will be only if the determinant of this system vanishes can be solved 
by diagonalizing of the matrix for E. The system ([I6D is called CI. This procedure goes 
far outside of the perturbation method. 

The CI approximation turned out to be extremely successful in atomic physics. In 
calculating atomic structure it allowed to increase the precision of calculating energy levels 
by one order. In the scattering processes it allowed to reveal fine resonance structure in 
scattering cross-sections due to formation of auto-ionizing states. So we expect that its 
applications will be even more important in strong interaction, where the perturbation 
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method is evidently not correct. The technique of apphcation of CIA is quite comphcated, 
since it needs to handle the matrices of large dimensions. In current work we used the 
code elaborated by O. Zatsarinny [|1^ . 

In this work we have applied the above described method for calculating P— wave 
"fine" -splitting of bb—, cc— and uu— systems to ^Pq; ^Pi and ^P2 levels. As well as for 
quarkonium constituted out of different masses like su, etc. In our case the corresponding 
operators His, Ht will be has the form 



LS 



1 1 

r 



12m2 



3— + {Ae - 1 ■ — e 

DVT 



S 



12- 



(19) 



It is important that all parameters except e are taken from [§,[1^, where excellent 
description of bottomonium and charmonium spectra were obtained. Moreover as it 
was shown in |jl3| the same parameters gives good p— meson trajectories. Actually the 
values ^ = 0.3 GeV^, p = 0.054 GeV were taken and as was taken in accordance to 
QCD. The only adjustable parameter was e. The experimental values are taken from 
||18|| . As mentioned above all calculations were carried out numerically. Special code was 
constructed for this purpose. The calculations were extended to fifth order in (plSj) ( see 
tables, rows 1-8), i.e. until the differences between the results did not go below several 
MeV level. 



Table 1. 66-system, 



0.3, £ = 0.5, TTib = 5.05 GeV. 



State 


AMth 


order 1 
MeV 


order 2 
MeV 


AMexp[18] 
MeV 


(0+ (2++))-X6i (0+(l++)) 


13P2 - l^Pi 


21.32 


21.23 


21.3 ± 1.3 


XM (0+(l++))-X6o(0+ (0++)) 


l^Pi - l3Po 


24.25 


25.51 


32.1 ±2 


Xt2 (0+(2++))-X6o(0+ (0++)) 


m - I'Po 


45.57 


46.74 


53.4 ± 1.9 


R 




0.88 


0.83 


0.66 


Xt2 (0+(2++))-X6i (0+(l++)) 


23p2 - 2^Pi 




17.72 


13.3 ±0.9 


Xn (0+(l++))-X6o(0+ (0++)) 


2-^Pi - 2^Po 




19.09 


23.1 ± 1.1 


Xb2 (0+ (2++))-X6o(0+ (0++)) 


23P2 - 23Po 




36.81 


36.4 ± 1.0 


R 






0.93 


0.57 


State 


AMth 


order 3 
MeV 


order 4 
MeV 


AMexpIIS] 
MeV 


X52 (0+ (2++))-X6i (0+(l++)) 


13P2 - l3Pi 


21.20 


21.19 


21.3 ± 1.3 


Xn (0+(l++))-X6o(0+ (0++)) 


l^Pi - l3Po 


26.08 


26.44 


32.1 ±2 


X,2 (0+(2++))-x.o(0+ (0++)) 


l'P2 - I'Po 


47.29 


47.63 


53.4 ± 1.9 


R 




0.81 


0.8 


0.66 


Xt2 (0+ (2++))-X6i (0+(l++)) 


23p2 - 2^Pi 


17.61 


17.58 


13.3 ±0.9 


XM (0+(l++))-x.o(0+ (0++)) 


2^Pi - 23Po 


20.37 


21 


23.1 ± 1.1 


Xt2 (0+ (2++))-X6o(0+ (0++)) 


23P2 - 23Po 


37.98 


38.58 


36.4 ± 1.0 


R 




0.86 


0.84 


0.57 
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Table 2. ^^^systcm, a, ^ 0.3, s = 0.45, nib = 5.05 GeV. 



State 


AMth 


order 1 
MeV 


order 2 
MeV 


AMexp[18] 
MeV 


X.2(0+(2++))-x.i(0+(l++)) 




20.00 


19.92 


21.3 ± 1.3 


(0+(l++))-X6o(0+ (0++)) 


V'P, V'Po 


23.37 


24.57 


32.1 ±2 


X;.2 (0+(2++))-x.o(0+ (0++)) 


1^P2 - l^Po 


43.37 


44.49 


53.4 ± 1.9 


R 




0.85 


0.81 


0.66 


Xb2 (0+ (2++))-X6i (0+(l++)) 


2^P2 - 2^ Pi 


— 


16.68 


13.3 ±0.9 


X61 (0+(l++))-X5o(0+ (0++)) 


2^ Pi - 23 Po 


— 


18.49 


23.1 ± 1.1 


(0+(2+-))-,uo(0+ (0+-)) 


2-^"P2 - 2--^Po 


— 


35.25 


36. 1± 1.0 


R 




— 


0.9 


0.57 


State 


AMth 


order 3 
MeV 


order 4 
MeV 


AMexpIIS] 
MeV 


X.2 (0+ (2++))-x.i (0+(l++)) 


13P2 - l^Pi 


19.89 


19.87 


21.3 lb 1.3 


X61 (0+(l++))-X5o(0+ (0++)) 


I'Pi 1-^Po 


25.13 


25.48 


32.1 ±2 


X/.2(0+(2+-))-AM-)(0+(0+-)) 


v"P2 r'^Po 


15.02 


15.36 


53. 1± 1.9 






0.79 


0.78 


0.66 


X.2 (0+ (2++))-Xw (0+(l++)) 


23p2 - 2^ Pi 


16.76 


16.65 


13.3 ±0.9 


X.i (0+(l++))-X6o(0+ (0++)) 


2^Pi - 2^Po 


19.73 


20.34 


23.1 ± 1.1 


X.2 (0+ (2++))-X(,o(0+ (0++)) 


2^P2 - 23Po 


36.38 


36.95 


36.4 ± 1.0 






0.85 


0.82 


0.57 



Table 3. cc-system, ag — 0.386, e = 0.3, rrie = 1.675 GeV. 



State 


AMth 


order 1 
MeV 


order 2 
MeV 


order 3 
MeV 


order 4 
MeV 


Xc2(0+(2++))-Xci(0+(l++)) 


13P2 - l3p, 


52.20 


51.66 


51.49 


51.41 


Xci (0+(l++))-Xco(0+ (0++)) 


l^Pi - l3Po 


69.31 


77.27 


81.55 


84.51 


Xc2(0+(2++))-Xc2(0+(0++)) 


1^P2 - I'PO 


121.51 


128.93 


133.04 


135.93 


R 




0.75 


0.66 


0.63 


0.61 


State 


AMth 


order 5 
MeV 


order 6 
MeV 


AMexp 
MeV 


[18] 




Xc2(0+(2++))-Xci(0+(l++)) 


V''P2 - V'Pl 


51.37 


51.33 


45.64 ±0.25 


Xci (0+(l++))-Xco(0+ (0++)) 


l^Pi - l^Po 


86.72 


88.48 


95.43 ± 1.12 


Xc2(0+(2++))-Xc2(0+(0++)) 




138.09 


139.82 


141.07 ± 1.13 


R 




0.59 


0.58 


0.48 
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Table 4. cc-system, a, ^ 0.386, e = OA, = 1.675 Gc^. 



State 


AMth 


order 1 
MeV 


order 2 
MeV 


AMexp[IS] 
MeV 


Xc2(0+(2++))-Xci(0+(l++)) 




66.81 


66.17 


45.64 ±0.25 


Xci (0+(l++))-Xco(0+ (0++)) 


r^Pi - r^Po 


79.36 


88.57 


95.43 ± 1.12 


Xc2(0+(2++))-Xc2(0+ (0++)) 




146.17 


154.74 


141. 07± 1.13 


R 




0.84 


0.75 


0.48 


State 




order 3 
MeV 


order 4 
MeV 


AMexp[1S] 
MeV 


Xc2(0+(2++))-Xci(0+(l++)) 


13P2 - l3Pi 


65.99 


65.92 


45.64 ± 0.25 


Xel (0+(l++))-Xco(0+ (0++)) 




93.45 


96.75 


95.43 ± 1.12 


Xc2(0+(2++))-Xc2(0+ (0++)) 


13P2 - l^Po 


159.45 


162.67 


141.07 ± 1.13 


R 




0.70 


0.68 


0.48 


Table 5. uu-system, — 0.52, e — 0.14, m„ = 0.33 GeV. 



State 


AMth 


order 1 
MeV 


order 2 
MeV 


AMexp[18] 
MeV 


a2(l-(2++))-ai(l-(l++)) 




15.34 


10.26 


88.2 ±40.7 


ai(l-(2++))-a„(l-(0++)) 


l^Pi - l3Po 


241.19 


316.91 


246.5 ± 40.9 


«2(l-(2+-))-«o(l-(0+-)) 


V''P2 I'Po 


256.53 


327.17 


33 1.6 ± 1.6 


R 




0.063 


0.032 


0.357 


State 


AMth 


order 3 
MeV 


AMexpIIS] 
MeV 




a2(l-(2++))-ai(l-(l++)) 


13P2 - l3Pi 


8.57 


88.2 ±40.7 


ai(l-(2++))-ao(l-(0++)) 


l^Pi - l3Po 


374.09 


246.5 ± 40.9 


a2(l-(2++))-ao(l-(0++)) 


1^P2 - l^Po 


382.66 


334.6 ± 1.6 


R 




0.023 


0.357 


Table 6. Mw-system, = 0.52, e = 0.145, m„ = 0.33 GeV. 




State 


AMth 


order 1 

MeV 


order 2 

MeV 


AMexp[18] 
MeV 


a2(l-(2++))-ai(l-(l++)) 


1^P2 - l^Pl 


24.41 


19.15 


88.2 ±40.7 


ai(l-(2++))-Oo(l-(0++)) 


l^Pi - l3Po 


247.52 


324.96 


246.5 ±40.9 


a2(l-(2++))-ao(l-(0++)) 


13P2 - l3Po 


271.92 


344.10 


334.6 ± 1.6 


R 




0.098 


0.059 


0.357 


State 


AMth 


order 3 

MeV 


AMexp[18] 
MeV 




a2(l-(2++))-ai(l-(l++)) 


1^P2 - l^Pl 


17.45 


88.2 ±40.7 


ai(l-(2++))-ao(l-(0++)) 


l^Pi - l3Po 


383.24 


246.5 ±40.9 


a2(l-(2++))-ao(l-(0++)) 


13P2 - l3Po 


400.69 


334.6 ± 1.6 


R 




0.046 


0.357 
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Table 7. gM-system, ct, = 0.42L e = 0.1875, m, = 0.5 GeV,mu = 0.33 GeV. 



State 


AMth 


order 1 
MeV 


order 2 
MeV 


order 3 
MeV 


order 4 
MeV 


AMexp[IS] 
MeV 


? 




43.16 


40.23 


39.29 


38.87 




? 




187.78 


230.62 


259.86 


282.63 




? 


l^P2 - I'Po 


230.94 


270.85 


299.15 


321.50 




R 




0.23 


0.17 


0.15 


0.14 





Table 8. g^-system, a, = 0.421, e = 0.2, = 0.5 GeV, ruu = 0.33 GeV. 



State 


AMth 


order 1 
MeV 


order 2 
MeV 


order 3 
MeV 


AMexp[18] 
MeV 


? 


13P2 - I'^Pi 


60.13 


56.95 


56.00 




? 


l^Pi - l^Po 


199.43 


244.92 


331.74 




■? 


r""P2 - i'''Pi> 


259.56 


301.81 


387.75 




R 




0.30 


0.23 


0.17 
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Let us make the following conclusions: 



1. The results for heavy quarkonium are quite good for values e = 0.3 — 0.45. Which is 



coincide with results obtained by Lai-Him Chan [|T9[ For light quarkonium the results are 
worse, which means that relativistic effects have to be taken into account more carefully. 
We believe, that the difference between e in case of bb and uu system exactly reflects this 
fact. 

2. The value of indicates that confinement has prevailingly scalar character. This 
conclusion does not contradict other authors pO|] . 

3. As it follows from (|18|) at e = 0.25 [Ae — 1) = the contribution of confinement 
vanishes totally. May be exactly this circumstances was the reason that some authors 
stated the pure one-gluon character of LS-splitting. 

4. The first column in tables corresponds to pure perturbation approach. It is clearly 
seen that this approach gives only rough qualitative estimate, but the results are dras- 
tically improving with switch on the CIA expansion. We believe that the using of CIA 
approach in quark physics will has bright future. Very convenient for fine-splitting char- 
acteristic is the coefficient 

M (3Pi) - M (3Po) ■ 

As cited in the experimental values of this parameter are R = 0.66 for bb, R = 0.48 
for cc and also R = 0.21 for uu. It is interesting to note that firstly this parameter was 
introduced in atomic physics where its value was established for atoms with two outer 
electrons (like Mg) as being R = 2. Since the fine mass splitting in atoms is totally due 
to one-photon exchange this value gives good indication as to the nature of the Lorentz 
character of qq— potential. As our calculations show the one-gluon term is totally vector 
while non-perturbative confinement term is of about 75-80% scalar. 

Confinement gives 20-25% plus to the results, obtained with using only the one gluon 
exchange term, for the description of the LS—, ST— interaction. Finally we obtain 
description of fine-effects at the quark systems with the precision 80-90%. 

These conclusions confirms with the qualitative estimation of Lucha and Shoberl . 



This result approximately confirms with the conclusion of Franzini |]I6[ in the part, that 
the confinement is pure scalar. We show that the results are very sensitive to the exact 
value of e and small change of e can destroy the agreement with the experiment what 
actually occurred in the icase of Franzini. 

We believe that exactly this feature will be successfully exploited in near future not 
only in hadronic structure, but also explaining resonances in vrvr— or tcN— resonances. 
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